**********************************************************
**********************************************************
** Replication File for Study 2
**
** Financial rewards do not stimulate co-production: 
**  Evidence from two experiments
**
** Voorberg, Jilke, Tummers, Bekkers
** 
** Public Administration Review
**
**********************************************************
**********************************************************


**Housekeeping
clear all
version 15
set more off


**Set working directory
cd "X" //Insert working directory for X


**Call data
use Study2-Data


**Produce Table 2
tab1 Gender Age educ Dutch,missing //educ: 1=lower, 2=higher, 3=uni, 4=other; Age: 1=<20, 2=20-39, 3=40-64, 4=65-80, 5=>80


**Mean differenece tests for DV1
ttest DV1 if Treat!=2,by(Treat) // Control vs Treat1
esize two DV1 if Treat!=2,by(Treat) // Cohen's d

ttest DV1 if Treat!=1,by(Treat) // Control vs Treat2
esize two DV1 if Treat!=1,by(Treat)

ttest DV1 if Treat!=0,by(Treat) // Treat1 vs Treat2
esize two DV1 if Treat!=0,by(Treat)

**Produce Figure 5
set scheme lean2
cibar DV1,over1(Treat) // Make sure to pre-install the "cibar" package


**Chi square tests for DV2
tab DV2 Treat if Treat!=2,col chi2 // Control vs Treat1
tab DV2 Treat if Treat!=1,col chi2 // Control vs Treat2
tab DV2 Treat if Treat!=0,col chi2 // Treat1 vs Treat2

**Produce Figure 6
cibar DV2,over1(Treat)



**Heterogenous treatment effects of DV1

**Gender
ttest DV1 if Treat!=1 & Gender==1,by(Treat) // Control vs Treat2 for Male only
esize two DV1 if Treat!=1 & Gender==1,by(Treat) // Cohen's d

ttest DV1 if Treat!=1 & Gender==2,by(Treat) // Control vs Treat2 for Female only
esize two DV1 if Treat!=1 & Gender==2,by(Treat)


**Age
ttest DV1 if Treat!=1 & Age==1,by(Treat) // Control vs Treat2 for <20 years
esize two DV1 if Treat!=1 & Age==1,by(Treat)

ttest DV1 if Treat!=1 & Age==2,by(Treat) // Control vs Treat2 for 20-39 years
esize two DV1 if Treat!=1 & Age==2,by(Treat)

ttest DV1 if Treat!=1 & Age==3,by(Treat) // Control vs Treat2 for 40-64 years
esize two DV1 if Treat!=1 & Age==3,by(Treat)

ttest DV1 if Treat!=1 & Age==4,by(Treat) // Control vs Treat2 for 65-80 years
esize two DV1 if Treat!=1 & Age==4,by(Treat)

ttest DV1 if Treat!=1 & Age==5,by(Treat) // Control vs Treat2 for 80+ years
esize two DV1 if Treat!=1 & Age==5,by(Treat)


**Place of residence
reg DV1 i.Treat##c.urban


**Confidence in government
reg DV1 i.Treat##c.trust


**Ethnicity
ttest DV1 if Treat!=1 & Dutch==0,by(Treat) // Control vs Treat2 for non-Dutch only
esize two DV1 if Treat!=1 & Dutch==0,by(Treat)

ttest DV1 if Treat!=1 & Dutch==1,by(Treat) // Control vs Treat2 for Dutch only
esize two DV1 if Treat!=1 & Dutch==0,by(Treat)


**Education
ttest DV1 if Treat!=1 & educ==1,by(Treat) // Control vs Treat2 for lower professional educ only
esize two DV1 if Treat!=1 & educ==1,by(Treat)

ttest DV1 if Treat!=1 & educ==2,by(Treat) // Control vs Treat2 for higher professional educ only
esize two DV1 if Treat!=1 & educ==2,by(Treat)

ttest DV1 if Treat!=1 & educ==3,by(Treat) // Control vs Treat2 for university educ only
esize two DV1 if Treat!=1 & educ==3,by(Treat)

ttest DV1 if Treat!=1 & educ==4,by(Treat) // Control vs Treat2 for other educ only
esize two DV1 if Treat!=1 & educ==4,by(Treat)


clear all
exit
